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Abstract 

We consider a system in which two users communicate with a destination with the help of a half- 
duplex relay. Based on the compute-and-forward scheme, we develop and evaluate the performance of 
coding strategies that are of network coding spirit. In this framework, instead of decoding the users' 
information messages, the destination decodes two integer-valued linear combinations that relate the 
transmitted codewords. Two decoding schemes are considered. In the first one, the relay computes 
one of the linear combinations and then forwards it to the destination. The destination computes the 
other linear combination based on the direct transmissions. In the second one, accounting for the side 
information available at the destination through the direct links, the relay compresses what it gets using 
Wyner-Ziv compression and conveys it to the destination. The destination then computes the two linear 
combinations, locally. For both coding schemes, we discuss the design criteria, and derive the allowed 
sum-rate. Next, we address the power allocation and the selection of the integer-valued coefficients to 
maximize the offered sum-rate; an iterative coordinate descent method is proposed. The analysis shows 
that our coding schemes can offer substantial gains over standard relaying techniques. The results are 
illustrated through some numerical examples. 
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I. Introduction 

Network coding was introduced by Ahlswede et al. in [T) for wired networks. It refers to each intermediate node sending 
out a function of the packets that it receives, an operation which is more general than simple routing O, O- In linear network 
coding, intermediate nodes compute and send out linear combinations over an appropriate finite field of the packets that 
they receive. In general, the function does not need to be linear. Although they are generally suboptimal for general wireline 
networks, linear network codes have been shown optimum for multicasting |4l [5 1. Moreover they have appreciable features, 
in particular simplicity (e.g., see (6), (7) are references therein). For these reasons, most of the research on network coding 
has focused on linear codes. 

The development of efficient network coding techniques for wireless networks is more involved than for wired network 
coding, essentially because of fading, interference and noise effects. For general wireless networks, the quantize-map-and- 
forward scheme of (8) and the more general noisy network coding scheme of J9) can be seen as interesting and efficient 
extensions for wireless settings of the original network coding principle. However, quantize-map-and-forward and noisy 
network coding are based on random coding arguments. For wireless networks, efficient linear network coding techniques 
make use of structured codes, and in particular lattices ifTol . Lattices play an important role in network coding for diverse 
network topologies, such as the two-way relay channel 1111 , lfl2l . the interference-relay channel 1131 and others. 

Recently, Nazer and Gastpar propose and analyse a scheme in which receivers decode finite-field linear combinations of 
transmitters' messages, instead of the messages themselves. The scheme is called "Compute-and-forward" (CoF) 1131 . and 
can be implemented with or without the presence of relay nodes. In this setup, a receiver that is given a sufficient number 
of linear combinations recovers the transmitted messages by solving a system of independent linear equations that relate the 
transmitted symbols. Critical in this scheme, however, is that the coefficients of the equations to decode must be integer- 
valued. This is necessitated by the fact that a combination of codewords should itself be a codeword so that it be decodable. 
Lattice codes have exactly this property, and are thus good candidates for implementing compute-and-forward. 

Compute-and-forward is a promising scheme for network coding in wireless networks. However, the problem of selecting 
the integer coefficients optimally, i.e., in a manner that allows to recover the sent codewords from the decoded equations and, 
at the same time, maximizes the transmission rate is not an easy task. As shown by Nazer and Gastpar (131 , the compute- 
and-forward optimally requires a match between the channel gains and the desired integer coefficients. However, in real 
communication scenarios, it is unlikely that the channels would produce gains that correspond to integer values. The selection 
of which integer combinations to decode is then a crucial task to be performed by the receivers. While it can be argued that 
linear combinations that are recovered at the same physical entity can always be chosen appropriately, i.e., in a way enabling 
system inversion to solve for the sent codewords, selecting these linear combinations in a distributed manner, i.e., at physically 
separate nodes, is less easy to achieve. By opposition to previous works, part of this paper focuses on this issue. 

In this work, we consider communication over a two-user multiaccess relay channel. In this model, two independent users 
communicate with a destination with the help of a common relay node, as shown in FigureQ] The relay is assumed to operate 
in half-duplex mode. 
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A. Contributions 

We establish two coding schemes for the multiaccess relay model that we study. The first coding scheme is based on 
compute-and-forward at the relay node. The relay uses what it receives from the transmitters during the first transmission 
period to compute a linear combination with integer coefficients of the users' codewords. It then sends this combination to the 
destination during the second transmission period. In addition to the linear combination that it gets from the relay's transmis- 
sion, the destination recovers the required second linear combination from what it gets directly from the transmitters, through 
the direct links. If the set of integer coefficients that are selected at the relay and the destination are chosen appropriately, 
the destination can solve for the transmitted codewords. We consider individual power constraints at the transmitters and the 
relay, and analyse the sum rate offered by this coding scheme 1141 lfT31 . 

In the second coding scheme both required linear integer combinations of the users' codewords are recovered locally at the 
destination. More specifically, the relay quantizes its output from the users' transmission during the first transmission period 
using Wyner-Ziv compression lfl6l . In doing so, it accounts for the output at the destination during this transmission period 
as available side information at the decoder. Then, the relay sends the lossy version of its output to the destination during 
the second transmission period. The destination determines the two required linear combinations, as follows. It utilizes an 
appropriate combination of the output from the users' transmission during the first period and of the compressed version of 
the relay's output during the second period; from this combination, two independent linear combinations relating the users' 
codewords are recovered. 

For the two coding schemes, we target the optimization of the transmitters and the relay powers, and of the integer 
coefficients of the linear combinations to maximize the achievable sum rate. These optimization problems are NP hard. For the 
two coding schemes, we develop an iterative approach that finds the appropriate power and integer coefficients alternatively. 
More specifically, we show that the problem of finding appropriate integer coefficients for a given set of powers has the same 
solution as an approximated mixed integer quadratic programming (MIQP) problem with quadratic constraints. Also, we show 
that the problem of finding the appropriate power policy at the transmitters and the relay for a given set of integer coefficients 
is a non-linear non-convex optimization problem. We formulate and solve this problem through geometric programming and 
a successive convex approximation approach ifTTl . 

Our analysis shows that both schemes outperform known strategies for this model that do not involve forms of network 
coding, such as those based on having the relay implement classic amplify-and-forward, decode-and-forward or compress- 
and-forward relaying schemes. We illustrate our results by means of some numerical examples. The analysis also shows the 
benefit obtained from allocating the powers and the integer coefficients appropriately. 

B. Outline and Notation 

An outline of the remainder of this paper is as follows. Section|n]describes in more details the communication model that 
we consider in this work. It also contains some preliminaries on lattices and known results from the literature for the setup 
under consideration where the relay uses standard techniques. In Section|ni] we describe our coding strategies and analyse 
the sum rates that are achievable using these strategies. Section|lV]is devoted to the optimization of the power values and the 
integer-valued coefficients for an objective function which is the sum-rate. Section|V]contains some numerical examples, and 
Section|Vl]concludes the paper. 
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We use the following notations throughout the paper. Lowercase boldface letters are used to denote vectors, e.g., x. Upper 
case boldface letters are used to denote matrices, e.g., X. Calligraphic letters designate alphabets, i.e., X. The cardinality of 
a set X is denoted by \X\. For matrices, we use the notation X s R mx ™, m, n e N, to mean that X is an m-by-n matrix, i.e., 
with m rows and n columns, and its elements are real-valued. Also, we use X T to designate the n-by-m matrix transpose of 
X; and, if m = n, det(X) to designate the determinant of X. We use I n to denote the n-by-n identity matrix; and to denote 
a matrix whose elements are all zeros (its size will be evident from the context). Similarly, for vectors, we write x e A™, e.g., 
A = R or A = Z, to mean that x is a column vector of size n, and its elements are in A. For a vector x e R™, ||x|| designates 
the norm of x in terms of Euclidean distance; and for a scalar x e R, \x\ stands for the absolute value of x, i.e., \x\ = x if 
x > and \x\ = —x if x < 0. For two vectors x 6 R™ and y s R™, the vector z = x o y e R™ denotes the Hadamard product 
of x and y, i.e., the vector whose ith element is the product of the ith elements of x and y, i.e., zi = (x o y) 2 = xiyi. Also, 
we use Var(x) to denote the power of x i.e. E[j|x|| 2 ]. The Gaussian distribution with mean fi and variance a 2 is denoted by 
Af(/J,, a 2 ). Finally, throughout the paper except where otherwise mentioned, logarithms are taken to base 2; and, for ifl, 
log + (a;) := max{log(a;), 0}. 

II. Preliminaries and System Model 

In this section, we first recall some basics on lattices, and then present the system model that we study and recall 
some known results from the literature, obtained through classic relaying, i.e., amplify-and-forward, decode-and-forward 
and compress-and-forward. The results given in Section lTl-Cl will be used later for comparison purposes in this paper. 

A. Preliminaries on Lattices 

Algebraically, an n-dimensional lattice A is a discrete additive subgroup of R". Thus, if Ai 6 A and A2 e A, then (Ai + 
A2) e A and (Ai - A2) e A. For an n-dimensional lattice A 6 R", there exists (at least one) matrix G € R nx ™ such that any 
lattice point A £ A can be written as an integer combination of the columns of G. The matrix G is called generator matrix of 
A, and satisfies 

A= {A = zG:zsZ n }. (1) 

A lattice quantizer Q\ : R™ -> A maps a point x e R™ to the nearest lattice point in Euclidean distance, i.e., 

Qa( x ) = ar g mm II x - All . (2) 
A e A 

The Voronoi region V(A) of A 6 A is the set of all points in R™ that are closer to A than to any other lattice point, i.e., 

V(A) = {xeR" : Q A (x) = A}. (3) 

The fundamental Voronoi region V of lattice A is the Voronoi region V(0), i.e., V = V(0). The modulo reduction with respect 
to A returns the quantization error, i.e., 

[x]modA = x-Q A (x) 6 V. (4) 

The second moment a\ quantifies per dimension the average power for a random variable that is uniformly distributed over 
V, i.e., 

(Ta = ; — r / l|x|| 2 dx (5) 

A nVol(V) Jv " " 
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where Vol( V) is the volume of V. The normalized second moment of A is defined as 

G(A) 



2 



(6) 



Vol(V) 2 /"' 

A lattice A is said to be nested into another lattice Afine if A c Afine, i-e., every point of A is also a point of Afine- We refer 
to A as the coarse lattice and to Afine as the fine lattice. Also, given two nested lattices A c Afine, the set of all the points of 
the fine lattice Afine that fall in the fundamental Voronoi region V of the coarse lattice A form a codebook 



C = Afine n V = {x = A mod A, A e Afine}- 



The rate of this codebook is 



i?=-log 2 (|C|). 

n 

Finally, the mod operation satisfies the following properties: 

(PI) [[x] modA + y] mod A = [x + y] mod A, Vx6l",ye 
(P2) [fc([x]modA)]modA= [fcxjmodA, V k e Z, x e K" 
(P3) 7([x] modA) = [7x] mod7A, V 7 e R, xel". 



(7) 



(8) 



(9) 



B. System Model 

We consider the communication system shown in Figure [TJ Two transmitters A and B communicate with a destination 
D with the help of a common relay R. Transmitter A, and Transmitter B want to transmit the messages W a £ Wa, and 
Wi, e Wft to the destination reliably, in 2n uses of the channel. At the end of the transmission, the destination guesses the pair 
of transmitted messages using its output. Let R a be the transmission rate of message W a and R}, be the transmission rate of 
message W^. We concentrate on the symmetric rate case, i.e., R a = R b = R, or equivalently, \W a \ = |W&| = 2 2nR . We measure 
the system performance in terms of the allowed achievable sum rate R sum = R a + Rb - 2i2. Also, we divide the transmission 
time into two transmission periods with each of length n channel uses. The relay operates in a half-duplex mode. 




Destination 



Fig. 1. Multiple-access channel with a half-duplex relay 
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During the first transmission period, Transmitter A encodes its message W a £ [l,2 2nfl ] into a codeword x a and sends 
it over the channel. Similarly, Transmitter B encodes its message W b £ [1, 2 2nR ] into a codeword and sends it over the 
channel. Let y r and be the signals received respectively at the relay and at the destination during this period. These signals 
are given by 

y r = h ar X a + h br X b + Z r 

Yd = h ad x a + h bd x b + z d , (10) 

where h ad and h bd are the channel gains on the links transmitters-to-destination, h ar and h br are the channel gains on the 
links transmitters-to-relay, and z r and z d are additive background noises at the relay and the destination. 

During the second transmission period, the relay sends a codeword x r to help both transmitters. During this period, the 
destination receives 

Yd = h rd x r + z d , (11) 

where h rd is the channel gain on the link relay-to-destination, and z d is additive background noise. 

Throughout, we assume that all channel gains are real-valued, fixed and known to all the nodes in the network; and the 
noises at the relay and the destination are independent among each others, and independently and identically distributed (i.i.d) 
Gaussian, with zero mean and variance N. Furthermore, we consider the following individual constraints on the transmitted 
power (per codeword), 

E[||x a || 2 ] = nfi'aP < nP a , E[\\x b \\' 2 ] =n/3 b P <nP b , E[||x r || 2 ] = n^P < nP r , (12) 

where P a > 0, P b > and P r > are some constraints imposed by the system; P > is given, and j3 a , f3 b and f3 r are some 
scalars that can be chosen to adjust the actual transmitted powers, and are such that < \/3 a \ < yPafP, < \f3 b \ < \JP b /P 
and < |/3 r | < \JP r /P. For convenience, we will sometimes use the shorthand vector notation = [h ad , h bd ] T , h r = 
[h ar , h br ] T £ R 2 and f3 = [/3 a , P b , f3 r ] T £ R 3 , and the shorthand matrix notation H = [h d ; ] £ R 2x2 . Also, we will 
find it useful to sometimes use the notation f3 s to denote the vector composed of the first two components of vector f3, i.e., 
/3 S = \_fia, /3 b ] T - the subscript "s" standing for "sources". Finally, the signal-to-noise ratio will be denoted as snr = P/N in 
the linear scale, and by SNR = 101og 10 (snr) in decibels in the logarithmic scale. 

C. Sum Rates Achievable Through Classic Relaying 

In this section, we review some known results from the literature for the model we study. These results will be used for 
comparisons in SectionlVl 

1 ) Amplify -and- Forward: The relay receives y r as given by dlOb during the first transmission period. It simply 
scales y r to the appropriate available power and sends it to the destination during the second transmission period. That is, the 
relay outputs x r = 7 y r , with 7 = \Z/3 2 snr/(l + snr \(3 S o h r || 2 ). The destination estimates the transmitted messages from its 
output vectors (y d , y d ). Using straightforward algebra, it can be shown ifTSIl that this results in the following achievable sum 
rate 
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where the vectors are given by = [h id , hi r h rd ^ / '(V 'l + 7 2 |^rd| 2 )] T f° r i = a,b, and the maximization is over /3. 
The achievable sum rate f\3[ does not require the two users to transmit at the same rate. Recall that, for a symmetric rate point 
to be achievable, both transmitters must be able to communicate their messages with at least that rate. Under the constraint of 
symmetric-rate, it can be shown rather straightforwardly lfT3l that the following sum rate is achievable with the relay operating 
on the amplify-and-forward mode, 

max - min j log (det (l 2 + /3aSnr(h a ha ))) , 

log (det (l 2 + /3 b 2 snr(h b hfT))) , ~ log (det (l 2 + ^(hahl) + /^snr(hi,h b r ))) |. (14) 

2 ) Decode-and-Forward: At the end of the first transmission period, the relay decodes the message pair (W a , W b ) 
and then, during the second transmission period, sends a codeword x r that is independent of x a and x b and carries both 
messages. The relay employs superposition coding and splits its power among the two messages. It can be shown easily that 
the resulting achievable sum rate is given by 1791 

Rfum = max i min | log + snr ||/3 S o h r || 2 ) .log (l + snr ||/3 S o h d || 2 ) + log(l + snr \h rd \ 2 f3 2 ^) |, (15) 

where the maximization is over (3. Under the constraint of symmetric-rate, it can be shown rather straightforwardly |[T3l that 
the following sum rate is achievable with the relay operating on the decode-and-forward mode, 

#°m = max imin |i?(h r ), R(h d ) + ~ log (l + snr|/i rd | 2 /? 2 ) |, (16) 

where 

R(hi) =min | log (l + sm\h ai \ 2 /3^j , log (l + sm\h bl \ 2 ^ , i log (l + snr \\/3 s o hj| 2 ) j. (17) 

3) Coinpress-and-Forward: At the end of the first transmission period, the relay quantizes the received y r using 
Wyner-Ziv compression [14], accounting for the available side information at the destination. It then sends an independent 
codeword x r that carries the compressed version of yy. The destination guesses the transmitted messages using its output from 
the direct transmission along with the lossy version of the output of the relay that is recovered during the second transmission 
period. It can be shown that the resulting achievable sum rate is given by |fl9l , 1201 , 

R*um = max \r CF , (18) 
4 



where 



, (l + snr||/3 s o h d || 2 ) (l + D/N + snr||/3 s ° h r | 2 ) - snr 2 ((/3 s o h r ) T (/3 s o h d )) \ 
R = max log ( ITTDjN) ' ^ 



the maximization is over (3 S and D > 0, where D is the distortion due to Wyner-Ziv compression, which is given by, 

TV 2 (l + snr||/3 s o h r || 2 ) 7V 2 (snr(/3 s o h r ) T ((3 s ° h d )) 2 
\h r d\ 2 Pr \h r d\ 2 Pr (1 + snr|/3 s ° h^P) 

Under the constraint of symmetric-rate, it can be shown rather easily that the following sum rate is achievable with the relay 
operating on the compress-and-forward mode, 



d cf 1 • I i it, |2 fl 2 snr|/i ar | 2 £ 2 \ / 2 a 2 &m\h br \ 2 ^ \ 1 C f 

R sum = max - mm < log I 1 + snr\h ad \ j3 a + — — — 1 , log 11 + snr\h bd \ /3 b + — — I , -R 



(21) 
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III. Network Coding Strategies 

In this section, we develop two coding strategies that are both based on the compute-and-forward strategy of fT-H . The 
two strategies differ essentially through the operations implemented by the relay. In the first strategy, the relay computes an 
appropriate linear combination that relates the transmitters' codewords and forwards it to the destination. The destination 
computes the other required linear combination from what it gets through the direct links. In the second strategy the relay 
sends a lossy version of its outputs to the destination, obtained through Wyner-Ziv compression 1161 - The destination then 
obtains the desired two linear combinations locally, by using the recovered output from the relay and the output obtained 
directly from the transmitters. 

A. Compute-and-Forward at the Relay 

The following proposition provides an achievable sum rate for the multiaccess relay model that we study. 

Proposition 1: For any set of channel vector h = [h ar , hb r , h ac i, /i&d, h r( f\ T e E 5 , the following sum rate is achievable 
for the multiaccess relay model that we study: 

where the maximization is over (3 and over the integer coefficients keZ 2 and teZ 2 such that |det(k, t)| > 1. 

As we already indicated, in the coding scheme that we use for the proof of Proposition Q] the relay first computes 
a linear combination with integer coefficients of the transmitters codewords and then forwards this combination to the 
destination during the second transmission period. The destination computes another linear combination that relates these 
codewords using its output from the direct transmissions. With an appropriate choice of the integer-valued coefficients of the 
combinations, the destination obtain two equations that can be solved for the transmitted codewords. 

Proof: In what follows, we first describe the transmission scheme, the encoding procedures at the transmitters and the 
relay, and the decoding procedures at the relay and the destination. Then, we analyze the allowed sum rate. 

Let A be an n-dimensional lattice of second moment a\ = P and normalized second moment G(A). Let V be its 
fundamental Voronoi region. Also, let Afine 2 A be chosen such that the codebook C = Afine n V be of cardinality 2 2nR 1211 . 
Also, let k = [k a , kb] 6 Z 2 and t = [t a , £f>] e ^ 2 be given such that |det(k, t)| = \k a tb - kf,t a \ > 1. 
The encoding and transmission scheme is as follows. 

Encoding: Let (W a , Wb) be the pair of messages to be transmitted. Let u a , and u r be some dither vectors that are drawn 
independently and uniformly over V and known by all nodes in the network. Since the codebook C is of size 2 2nR = |W a |, 
there exists a one-to-one mapping function 4> a (-) between the set of messages {W a } and the nested lattice code C. Similarly, 
there exists a one-to-one mapping function 4>b{') between the set of messages {Wb} and the nested lattice code C. Let 
v a = 4>a(W a ) and v b = (f>b(W b ), where v a e C and v 6 e C. 

During the first transmission period, to transmit message Wa, Transmitter A sends 

Xa = Pa ([v a - u a ] mod A) , (23) 
for some /J a et such that < \J3 a \ < \fPajP\ and to transmit message Wb, Transmitter B sends 

x 6 = Pb ( [v fc - u 6 ] mod A) , (24) 
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where < \j3 b \ < \fPtJP- The scalars /3 a and j3 b are chosen so as to adjust the transmitters' powers during this period. 

As will be shown shortly the relay decodes correctly an integer combination e2 = k a v a + k b v b from what it receives during 

the first transmission period. It then sends 

x r = f5 r ( [k a v a + k b v b - u r ] mod A) (25) 

during the second transmission period, where the scalar /3 r is chosen so as to adjust its transmitted power during this period. 
Decoding: During the first transmission period, the relay receives 

y r = harX-a + h br x b + z r . (26) 
Next, the relay performs the following modulo-reduction operation on the received signal: 

y'r = [a r y r + k a u a + k b u b ] mod A 

= [av(/i or x a + h br x b + z r ) + ^x a - ^x„ + ■^■x b - -^-Xf, + k a u a + k b u b ] mod A 

Pa Pa Pb Pb 

^= [(a r h ar - 7r") x a + (a r h br - ■^■)x b + fc a [v Q - u a ] mod A + k b [v b - u b ] mod A + a r z r + k a u a + k b u b ] mod A 

Pa Pb 

( =' [k a v a + k b v b + (a r h ar - -^-)x a + {a r h br - ■^)x h + a r z r ] mod A 

Pa Pb 



(c) 



[fc a v a + fcfcVb + z' r ] mod A (27) 



where (a) follows by substituting x a and x b using i23l and (1241) : (6) follows by using the properties (PI) and (P2) in |[9); and 
(c) follows by substituting, 

i k a k b 

z'r = [a r z r + (a r h ar - -^)x a + (ce r h br - — )x h ] mod A. (28) 

Pa Pb 

The parameter a r e R is some inflation factor whose optimal value will be specified below and z' r is the effective noise at the 
relay. Since an integer combination of lattice points is a lattice point, e2 = [k a v a + k b v b ] s A and, thus, the equivalent channel 
y' r is a modulo-lattice additive noise (MLAN) channel 1221 with noise equal to z' r . By decoding the lattice point e2 e A from 
the MLAN channel y' r , the relay obtains an integer combination of the users' codewords, e-2 = [&aV a + k b v b ]. 
During the first transmission period, the destination receives 

Yd = KdX-a + h bd x b + z d . (29) 

Similar to the relay, the destination computes a linear combination with integer coefficients of the transmitters' codewords by 
performing the modulo-reduction operation: 

y'd = [adYd + tau a + t b u b ] mod A 

= [ta\a + t b V b + z' d ] mod A (30) 

where a d £ K is some inflation factor and z' d is the effective noise at the destination, given by 

z' d = [a d z d + (a d h ad - ^-)x„ + (a d h bd - ^-)x 6 ] mod A. (31) 

Pa Pb 

Thus, during the first transmission period, by using the MLAN channel y' d , the destination obtains a first integer combination 
ei = [t a v a + t b v b ] s A of the users' codewords. 
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During the second transmission period, the destination receives from the relay the signal, 

Yd = Kd*r + Z d 

= KdPr ([fc Q V a + k b V b - U r ] mod A) + Z d . (32) 

Again, by performing a modulo-reduction operation on the obtained signal, the destination gets 

Yd = [&dYd + Ur] mod A 

= [fcaVa + k b v b + z d ] mod A, (33) 

where a d £ R is some inflation factor whose optimal value will be specified below and z' d is the effective noise at the 
destination, given by 

Zd = [&did + (a d h r d ~ tt)^] mod A. (34) 

Pr 

Thus, during the second transmission period, the destination obtains a second linear combination e2 = [fc a v a + t b v b ] £ A that 
relates the users' codewords. This combination has been computed and forwarded by the relay. 

Summarizing, over the entire transmission time, the destination collects two linear combinations with integer coefficients 
that relate the users' codewords, as 



ei 




ta 


tb 




V a 


e2 






h 







(35) 



Now, since the integer-valued matrix in i35i is invertible (recall that the integer-valued coefficients are chosen such that 
det(k, t) + 0), the destination obtains the transmitted codewords by solving d35| i. 

The destination is able to recover the messages W a and W b reliably only if the message rate is less than or equal to the 
computational rate 1131 . 

Rate Analysis: 

Step 1) First, we find the rate at which the relay decodes reliably the combination e2 = [fc a v a + t b v b ] 6 A. Consider the 
channel from e2 to y' r as given by d27l >. Due to the dither u a , the input x a of Transmitter A is independent of v a = 4>a(W a ) £ 
C, and is uniformly distributed over V (see, e.g., l23l . 1241 ). Similarly, due to the dither u b , the input Xf, of Transmitter B is 
independent of v;, = 4> b (W b ) £ C, and is uniformly distributed over V. The effective noise z' r is then independent of v a and 
v , and so of e2 = fc a v a + k b v b . It is composed of a "self noise" component and a Gaussian noise component; and its average 
power satisfies 

Var (z r ) < Var(a r z r + (a r h ar - -^)x a + (a r h br - — )x b ) 

Pa Pb 
= a r N + P(a r harf3a - k a ) 2 + P{a r h br p b - k b ) 2 . (36) 

where the inequality follows since the modulo-reduction operation can only reduce the entropy. 
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The rate (per channel use) at which the relay decodes the linear combination with integer coefficients e 2 = k a v a + k b v b is 
given by (25), 

^(e 2 ;yr) = ^ h (y'r) - j-A(yrlea) 

( =^(y;)-f K<) 

2n 2n 

(e J -Llog + Vo\(V)-±h(z' r ) 

2n 2n 

( > - log + | ^5 — n -| - - log (27reG(A)) 

% ] I i og+ f fc 2 + k l _ ^y^/^y' - i log (2«G(A)) 



= 1 log + f ( ||kf - %^f| - - log(2. eG (A)) (37) 

where (d) follows from the fact that x' r is independent of e 2 ; (e) follows from the fact that y' r is uniformly distributed over V; 
(/) follows from the fact that the entropy of z' r is upper-bounded by that of an i.i.d Gaussian vector that has the same second 
moments and by using J36l ); and (g) follows through straightforward algebra and by noticing that the right hand side (RHS) 
of d36b is maximized with the value of a r set to 

^° h ") Tk (38 ) 



N + P\\/3 3 oh r || 2 ' 

Step 2) Using analysis and algebra that are essentially similar to those in Step 1) above, it can be shown that the destination 
decodes the linear combination ei = t a v a + t b v b from what it gets directly from the direct links at rate 

i w<> , \ * ((w - 7f;;:?S )~') - 3 '°*< 2 *« g(a » ■ 

Similarly, it can be easily shown that the destination decodes the linear combination that is computed at the relay and forwarded 
to the destination during the second transmission period at rate 

^/(e 2 ; y' d ) > \ log (l + PlK f p2r ) - \ log (27reG(A)) . (40) 

Let us define R S1 (A) as the RHS of 07), 7? sd (A) as the RHS of lO, and 7? rd (A) as the RHS of J40t . 

The above means that using the coding scheme that we described, which employs the n-dimensional lattice A, the 
destination can decode the transmitters' codewords correctly at a transmission rate equal to the minimum among i? s ,(A), 
-Rsd(A) and -R rd (A), i.e., R(A) = min{i? sl (A), i? sd (A), i? rd (A)}. The allowed sum rate is given by i^(A) = 2i?(A) = 
2 min{i? sr (A), _R sd (A), _R rd (A)}. Noticing that 27reG(A) -> 1 when n -> oo |1Q| . the desired sum rate l l22l ) is obtained by 
taking the limit of R^°m (A) as n goes to infinity; and this completes the proof of Proposition Q] □ 

Although the achievable sum rate in Proposition[T]requires the relay to only decode a linear combination of the codewords 
transmitted by the users, not the individual messages, this can be rather a severe constraint in certain cases. In the following 
section, the relay only compresses its output and sends it to the destination. The computation of the desired linear combinations 
of the users' codewords takes place at the destination, locally. 
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B. Compress-and-Forward at the Relay and Compute at the Destination 

The following proposition provides an achievable sum rate for the multiaccess relay model that we study. 

Proposition 2: For any set of channel vector h = [h ar , hb r , h ac [, hf, d , h rd ] T e R 5 , the following sum rate is achievable: 

„CoD 1 ■ Jl +/ Sm " ^ 

R sum = max - mm \ log — — = — — — , 

2 \ snr /3 s oH i a t -t 2 + (a t oa t ) T n d ) 



lo +/ snr \ 

\ snr||/3 s o H T a k - k|| 2 + (a k o a k ) T n d ) 



(41) 



where at = [ait, &2t\ an d a k = [ai k , a 2k ] £ K are some inflation factors, = [1, 1 + D/N] e R , and D is given by 

d _ 7V 2 (l + snr||/3 s oh r | 2 ) 7V 2 (snr(/3 s o h r ) r (/3 s o h d )f 
\h rd \ 2 P r \h r d\ 2 Pr (1 + snr||/3 s o h d | 2 ) 

and the maximization is over at, a k , f3 s , and over the integer coefficients k and t such that |det(k, t)| > 1. 

As we indicated previously, in the coding scheme that we use for the proof of Proposition [2] the relay conveys a lossy 
version of its output to the destination during the second transmission period. In doing so, it accounts for the available side 
information at the destination, i.e, what the destination has received during the first transmission period. The destination 
computes two linearly independent combinations that relate the users' codewords using its outputs from both transmission 
periods, as follows. The destination combines appropriately the obtained lossy version of the relay's output (that it recovered 
from the relay's transmission during the second transmission period) and from what it received during the first transmission 
period. Then it computes two linearly independent combinations with integer coefficients that relates the users' codewords. 

Proof: The transmission scheme, and the encoding procedures at the transmitters are similar to those in the proof of 
Proposition Q] Therefore they will be outlined only, for brevity. We will insist more on aspects of the coding scheme that are 
inherently different from those of the coding scheme of Proposition Q] 

Encoding: During the first transmission period, the transmitters send the same inputs as in the coding scheme of Proposi- 
tion!]] i- e -> to transmit message W a , Transmitter A sends the input x a given by 1 123k and to transmit message Wt,, Transmitter 
B sends the input X£, given by ll24l . 

The relay quantizes what it receives during the first transmission period using Wyner-Ziv compression lfl6l . accounting for 
the available side information y d at the destination. Let y> be the compressed version of yy given by 

yy = y r + d (43) 

where d is a Gaussian random vector whose elements are i.i.d with zero mean and variance D; and is independent of all 
other signals. Also let _Rwz be the resulting compression rate. During the second transmission period, the relay conveys the 
description y r of yy to the destination. To this end, it sends an independent Gaussian input x r with /? 2 P and carries the 
Wyner-Ziv compression index of yy. 

Decoding at the destination: During the two transmission periods, the destination receives, 

Yd = h ad x a + h bd x b + z d 

Yd = h rd x r + z d . (44) 
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The destination computes two linearly independent combinations with integer coefficients that relate the users' codewords, 
as follows. 

Step 7) It first recovers the compressed version of the relay's output sent by the relay during the second transmission period, 
by utilizing its output y d as well as the available side information y d . As it will be shown below, the destination recovers the 
compressed version y r of y r if the constraint l f54l below is satisfied (see the "Rate Analysis" section). 
Next, the destination combines y d and y r and uses the obtained signal to compute a linear combination with integer coeffi- 
cients of the users' codewords. More specifically, let 

yt = oL lt y d + £*2tyy 

= \a\thad + a 2t h a r)x-a + {ct\th bd + a 2t h br )ic b + a u z d + ot, 2t z r + a 2t d, (45) 

for some at = [ait,a 2 t] T £ K 2 - The destination uses the obtained signal yt to compute a linear combination with integer 
coefficients of the transmitters' codewords by performing the modulo reduction operation 

Yt = [yt + taVLa + t b u b ] mod A 

= [£ a v a + t b v b + z' t ] mod A (46) 

where the algebra follows l l27l > and z' t is the effective noise given by 

z' t = \a u z d + a 2t z r + a 2t d + (ai t h ad + a 2 th a r - + {a\ t h bd + a 2f h br - -^-)x b l mod A. (47) 

L Pa Pb J 

Finally, by decoding the lattice point ei = [t a v a + t b v b ] e A using the MLAN channel y' t , the destination obtains a first linear 

combination with integer coefficients of the users' codewords. 

Step 2 Second, the destination again combines y^ and y r and uses the obtained signal to compute a second linear 

combination with integer coefficients of the users' codewords which is different from the one decoded in Step 1. More 

specifically, let 

Yfc = aikYd + a 2 kfr 

= (&lkhad + a 2k h ar )y^a + {aikHd + a 2k h br )yi b + a lk z d + a 2k z r + a 2k d, (48) 

for some a k = [ai k ,a 2k ] T e R 2 . The destination uses the obtained signal y^ to compute a linear combination with integer 
coefficients of the transmitters' codewords by performing the modulo reduction operation 

y'k = [yk + k a u a + k b u b ] mod A 

= [fc a v a + k b v b + z' k ] mod A (49) 

where the algebra follows l l27l > and z' k is the effective noise given by 

Zfc - \ai k z d + a 2k z r + a 2k d + (ai k h ad + a 2k h ar - 7r)x a + i^lkHd + a 2k h br - 7p)xJ mod A. (50) 

L Pa Pb 1 

Finally, by decoding the lattice point e 2 = [fc a v a + k b v b ] e A using the MLAN channel y' k , the destination obtains a second 
linear combination with integer coefficients of the users' codewords. 
Rate Analysis: 
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The relay compresses its output y r at the per-channel use Wyner-Ziv quantization rate 

Pwz = r--f(yr;yr|yd) 

2n 

= -^-KyAyd) - -^-h(fr\yr,yd) 

2n 2n 

2n 2n 

{b = ] -^log(27r e ) n |E[ yi .y^] -E[y r E[y r | y(i ] T ]| - \ log (JhreD) 
4n I 14 

(c)l / iV + P|/3 s °h r | 2 [P(/3 s oh r ) T (/3 s oh rf )] 2 ^ 

= I log l 1 + 5 ^vTpP^P) J' (51) 

where (a) follows since d is independent of (y?-,yd); (b) follows since (yv,yd) are jointly Gaussian and d is a Gaussian 
vector; (c) follows through by straightforward algebra, and by noticing that the minimum mean square error (MMSE) 
estimator of y r given y d is given by 

F r. I 1 P(/3 s oh r ) T (/3 s oh d ) 

E[yrlyd] = nTpW^p yd - (52) 

At the end of the second transmission period, the destination can decode the correct relay input x r reliably if 

2 a 2 • 



From < I5 1 b and | |53"1 ). we get the following constraint on the distortion 



Pwz < -log|l+ 1 7; rr |. (53) 



N 2 (l + snr||/3 s o h r || 2 ) 7V 2 (snr(/3 s o h r ) T ((3 s o h d )) 2 
D> — — — — — — . (54) 

\h r d\ Pr \h r d\ 2 Pr(l + snr||/3 s o h d || 2 ) 



The above implies that, under the constraint I I541 ). the destination recovers the lossy version y r of what it was sent by the relay 
during the second transmission period. 

The destination processes yt and y^ to obtain the linear combinations ei = [t a v a + tfev^] , and e2 = [fc a v a + k^vi,] of 
the users' codewords. Using the MLAN channel y' t given by l B6l > and proceeding similar to above, we get that the destination 
obtains the desired combination t at mutual information satisfying 

^-/(ei;y£)> -\og + ( ^— )_ii g( 27re G(A)), (55) 

2n 4 \sm\\/3 s 0H 1 at - t\\ z + (at o at) 1 n d / 4 

where the distortion D satisfies the constraint fl54b and at should be chosen to minimize the effective noise z' t in BTl ). i.e., 

such that 

at = (GG T + N d ) _1 Gt, (56) 

where G = [(0 S o h d ) T ; (/3 S o h r ) T ] e R 2x2 and N d = [1/snr, 0; 0, 1/snr + D/P] e R 2x2 . In a similar way, we get that the 
destination obtains the desired combination k at mutual information satisfying 

^/(e 2 ;yfc)> ]log + ( )_Ii og ( 27 r e G(A)), (57) 

2n 4 \snr||/3, o H J a fe - k|| 2 + (a fc o a fc ) J n d / 4 

where the distortion D satisfies the constraint d54| i, and aj. should be chosen to minimize the effective noise z' k in d50l >, i.e., 

such that 

a* k = (GG T + N d ) _1 Gk. (58) 
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Let us define J?i(A) as the RHS of (|55) and R 2 (A) as the RHS of {57). Tne above 

means that using the lattice-based 

coding scheme that we described, the destination can decode the transmitters' codewords correctly at the transmission sum 
rate R^"m (A) = 2 min{i?i(A), i?2(A)} provided that the condition ([54) is satisfied. Furthermore, investigating the expression 
of i?i(A), it can easily be seen that it decreases with increasing D. Also, observing that the RHS of d54l decreases if /3 r 
increases, the largest rate -Ri(A) is then obtained by taking the equality in the distortion constraint d54t with /3 r = P r /P- 
Finally, observing that 2-7reG(A) -> 1 when n -> oo 1101 , the desired sum rate i ]41 1 is obtained by taking the limit of Rf^ (A) 
as n goes to infinity; and this completes the proof of Proposition [2] □ 
Remark 1: There are some high level similarities among the coding strategies of proposition [TJ and proposition |7J In 
particular, they both consist essentially in decoding two linearly independent equations. However, as we mentioned previously, 
the required two equations are obtained differently in the two cases. More specifically, while the two equations are computed 
in a distributed manner using the coding strategy of proposition [TJ they are both computed locally at the destination in a 
joint manner using the coding strategy of proposition [2] A direct consequence of all the processing being performed locally 
at the destination with the latter coding scheme is that both computations of the required equations utilize all the output 
available at the destination, i.e., the output received during the first transmission period as well as the output received during 
the second transmission period, in a joint manner. Recall that, by opposition, the coding strategy of propositionQ]is such that 
the computation of one equation utilizes only the output received directly from the transmitters during the first transmission 
period, and the computation of the other equation is limited by the weaker output among the output at the relay during the 
first transmission period and the output at the destination during the second transmission period (since the equation decoded 
at the relay has to be recovered at the destination). The joint processing at the destination somehow gives some advantage 
to the coding strategy of proposition [2] The reader may refer to Section [V] where this aspect will be illustrated through some 
numerical examples and discussed further. 

IV. Sum Rates Optimization 

Section IIV-AI is devoted to find optimal powers and integer-coefficients that maximize the sum rate of Proposition [TJ 
Section ITV-B I deals with the optimization problem of Proposition [2] 

A. Compute-and-Forward at Relay 

1 ) Problem Formulation: Consider the sum rate Rf°m as given by ((22) in Proposition!]] The optimization problem 
can be stated as: 

(A) : max | min {log- ((||tf - Xlff h J , ^ ((||k|| 2 - WpftStf f) ■ (* + ^)} <«» 

where the maximization is over (3 such that < \/3 a \ < yPa/P, < |/3f,| < \fPbfP, < |/3 r | < \J P r /P and over the integer 
coefficients k and t such that |det(k, t)| > 1. 

The optimization problem (A) is non-linear and non-convex. Also, it is a MIQP optimization problem; and, so, it is not 
easy to solve it optimally. In what follows, we solve this optimization problem iteratively, by finding appropriate preprocessing 
vector (3 and integer coefficients t and k alternatively. We note that the allocation of the vector j3 determines the power that 
each of the transmitters and the relay should use for the transmission. For this reason, we will sometimes refer loosely to the 
process of selecting the vector f3 appropriately as the power allocation process. 
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Let, with a slight abuse of notation, denote by i? su ° m [t] the value of the sum rate at some iteration i > 0. To compute R s °m 
as given by l |59t iteratively, we develop the following algorithm, to which we refer to as "Algorithm A" in reference to the 
optimization problem (A). 



Algorithm A Iterative algorithm for computing B^°^ as given by d59b 
1: Initialization: set i = 1 and (3 = . 

2: Set f3 = (3^~ 1 ' in ( 1591 ), and solve the obtained problem using Algorithm A-l given below. Denote by k^'' the found k, and 
by t (,) the found t. 

3: Set k = k*- 1 ' and t = t^ 4 ' in and solve the obtained problem using Algorithm A-2 given below. Denote by /3^' the 
found (3. 

4: Increment the iteration index as i = t + 1, and go back to Step 2. 

5: Terminate if ||/9« - ^ \\ < e, \R^[l] - R^[l - 1]| < e or i&] < R^[l - 1]. 



As described in Algorithm A, we compute the appropriate preprocessing vector f3 and integer coefficients k and t, 
alternatively. More specifically, at iteration i > 1, the algorithm computes appropriate integer coefficients k^ e 1? and 
t/ 1 - 1 s Z 2 that correspond to a maximum of ( |59t computed with the choice of the preprocessing vector (3 set to its value 
obtained from the previous iteration, i.e., (3 = /3^ t ~ 1 ' ) (for the initialization, set (3^ to a default value). As we will show, this 
sub-problem is a MIQP problem with quadratic constraints; and we solve it iteratively using Algorithm A-l. Next, for the 
found integer coefficients, the algorithm computes adequate preprocessing vector f3^ that corresponds to a maximum of l |59t 
computed with the choice k = and t = t' l \ As we will show, this sub-problem can be formulated as a Complementary 
geometric programming problem. We solve it through a geometric programming and successive convex optimization approach 
(see Algorithm A-2 below). The iterative process in Algorithm A terminates if one of the following two conditions holds: 1) 
either of j/3' 1 ' - f3^ l ~ 1 ^ || and \R^m [*■] _ ^Snt 1 — 1]| is smaller than a prescribed small strictly positive constant e — in this 
case, the optimized value of the sum rate is i?^ u °m[t], and is attained using the preprocessing power vector (3* = (3^' and 
integer vectors k* = k^ and t* = t^ u '\ or 2) the sum rate stops increasing, i.e., ^ -Rsum[> ~ 1] — m this case, the 

optimized value of the sum rate is -R£°m [t - 1], and is attained using the preprocessing power vector /3* = f3^ L ~ 1 ' and integer 
vectors k* = k (,_1) and t* = t (t_1) . 

In the following two sections, we study the aforementioned two sub-problems of problem (A), and describe the algorithms 
that we propose to solve them. 

2) Integer Coefficients Optimization: In this section, we focus on the problem of finding appropriate integer 
vectors keZ 2 and teZ 2 for a given choice of the preprocessing vector (3. Investigating the objective function in (3§), it can 
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easily be seen that this problem can be equivalently stated as 

Ai 



mm 
k,t, A t 



s. t. 



Ai > iitir ■ 



Ai > Ilk 



P((/3 s °h d ) T t) 2 
iV + P||/3 s °h d p 

2 P((/3 s oh r ) T k) 2 

iV + P||/3 s oh r || 2 
N 



Al > N + $\h rd \ 2 P' 
|det(k,t)| > 1 



keZ 2 , tsZ 2 , Ai s: 



(60a) 
(60b) 

(60c) 

(60d) 
(60e) 
(60f) 



Note that Ai is simultaneously an extra optimization variable and the objective function in ( 1601 ). Also, it is easy to see that 
the integer coefficients k and t that achieve the minimum value of Ai also achieve a maximum value of the objective function 
in (HI. 

In order to reformulate problem d60b in a manner that will be convenient for solving it, we introduce the following quanti- 
ties. Leta = [0,0,0,0, 1] T ; ai = a 2 = a 3 = [0, 0, 0, 0, -1] T and a 4 = [0, 0, 0, 0, 0] T . Also, let b = [t a , t b , k a , k b , A X ] T ; and 
the scalars c\ = c 2 = 0, C3 = N/(N + /3r\h rdi \ 2 P) and C4 = -1. We also introduce the following five-by-five matrices Fi, F2, 
F3 and F4, where 

r 

2(i 2 -n 2 ) 



Fi 



2(I 2 -f2i) 




F 2 















(61) 



with 



Sli := 

n 2 := 



N + P\\h d \\ 2 
P 

iV + P||h r || 2 



(p s oh d )(i3 s °h d y 



(/3 s oh r )(/3 s oh r ) J , 



0-20 
\ = 2 


The optimization problem ( I60t can now be reformulated equivalently as 



(62) 



(63) 



mm 
b 



s. t. 



a b 

It t 

-b F 2 ;b + &i b < Cj i = l,—,4 
k 6 Z 2 , t £ Z 2 , Ai 6 R 



(64) 



The equivalent optimization problem i64l is a MIQP problem with quadratic constraints (26). If the involved matrices 
associated with the quadratic constraints (i.e., the matrices Fi, F 2 and F4 here) are all semi-definite, there are known 
approaches for solving MIQP optimization problems, such as cutting plane, decomposition, logic-based and branch and bound 
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approaches (26). In our case, it is easy to see that the matrices Fi and F2 are positive semi-definite. However, the matrix F4 
is indefinite, irrespective to the values of k and t. 

In order to transform the optimization problem $64[ into one that is MIQP-compatible (i.e., in which all the quadratic 
constraints are associated with semi-definite matrices), we replace the quadratic constraint J60et with one that is linear, as 
follows. We introduce the following two real-valued vectors k = [k a , k b ] T e R 2 and t = [t a ,t b \ T £ R 2 defined such that they 
satisfy 

k = /«oexp(k), t = roexp(t), (65) 



where k = [k q , K b ] £R and r = [r a , r b ] £R are constant vectors to be chosen appropriately. Thus, the constraint d60eb 
can now be rewritten equivalently as 

|det(k,t)| > 1 iff \n a T b exp(k a + i b ) - K b T a exp(k b + t a )\ > 1. (66) 

Now, we linearize the constraint ( 1661 ) by selecting the constant vectors n and r such that the first order Taylor series 
approximations exp(k) rj 1 + k and exp(t) « 1 + t hold. Hence, the constraint i |60e| > can be rewritten as 



|det(k,t)|>l iff \n a T b (l + k a +t b ) - K, b T a (l + k b + t a )\>l. (67) 

Note that the constraint fl67| > is now linear. The optimization problem d64l > has the same solution as the following problem 
which is MIQP-compatible, 

min Ai (68a) 

k,t, Ai 

s.t. itf _^;M T t) 2 _ Ai<0 

2 _p ( (/3 oh r rk) 2 

" " iV + Pl^olvl 2 

N 

Ai < (68d) 



N + $\h rd \ 2 P 

-|K Q r fc (l + k a + i b ) - n b T a (l + k b + t a )\ < -1 (68e) 
k- k ~ 

— - 1 - hi < 0, — - + 1 + ki < 0, i = a, b (68f) 

— - 1 - ii < 0, -— +l+k<0, i = a,b (68g) 

n n 

k, t s Z 2 , k, t, k, t e R 2 , Ai 6 R. (68h) 

The optimization problem l l68l can be solved iteratively using Algorithm A-l hereinafter. 

3) Power Allocation Policy: Let us now focus on the problem of finding an appropriate preprocessing vector (3 
for given integer vectors k £ Z 2 and t e Z 2 . Again, investigating the objective function in d59t . it can easily be seen that this 
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Algorithm A-l Integer coefficients selection for as given by ( |59l l 
1: Initialization: set ti = 1. 

2: Use the branch-and-bound algorithm of 1271 . 1281 to solve for Aj k*-' 1 ' and t < - L1 ' > with the constraint d68et substituted 

With -K a T b (l + k a + tb) + K b T a (l + k b + t a ) < -1. 

3: Update the values of k and t in a way to satisfy d65| l; and increment the iteration index as i\ = ti + 1. 

4: Terminate if exp(k^ 11 ^) 1 + jo'" and exp(t^ 11 ') « 1 H-t' 11 '. Denote the found solution as A™" 1 ' 1 . 

5: Redo steps 1 to 4 with in Step 2 the constraint ( |68et substituted with K a T b (l + k a + t b ) - K b T a (l + k b + t a ) < -1. In Step 

4, denote the found solution as A™"' 2 . 
6: Select the integer coefficients corresponding to the minimum among A™"' 1 and A™"' 2 . 



problem can be equivalently stated as 



mm 

P, A 2 



s. t. 



A 2 

A 2 > ||t| 

A 2 > ||k 
A 2 > 



a P((/3 s oh d ) T t) 2 
Ar + p||/3 s oh d p 

2 P((/3 s oh r ) T k) 2 

" iV + P||/3 s oh r || 2 
N 



N + P\h rd \ 2 l3?' 



<Pi< 



f3 s eR 2 , /3sK 3 , A 2 e: 



(69a) 

(69b) 

(69c) 
(69d) 

(69e) 
(69f) 



Here, similarly to the previous section, A 2 is simultaneously an extra optimization variable and the objective function in 
d69b . Also, it is easy to see that the value of (3 that achieves the minimum value of A 2 also achieves a maximum value of the 
objective function in l l59t . 

The optimization problem in l |69l is non-linear and non-convex. We use geometric programming 1 1 71 for solving it. Geometric 
programming is a special form of convex optimization for which efficient algorithms have been developed and are known in 
the related literature ||29| . There are two forms of GP: the standard form and the convex form. In its standard form, a GP 
optimization problem is generally written as ||29l 



minimize 



subject to 



/o(/3,A 2 ) 
/,(/3,A 2 )<l, 
ffz(/3,A 2 ) = l, 



J = l 
1 = 1, 



,L, 



(70a) 
(70b) 
(70c) 



where the functions fa and fj, j = 1, J, are posynomials and the functions gi, I = 1, L, are monomials in j3 and A 2 . In 
its standard form, d70fc is not a convex optimization problem. However, when possible, a careful application of an appropriate 
logarithmic transformation of the involved variables and constants generally turns the problem d70t into one that is equivalent 
and convex. That is, {70} is a GP nonlinear, nonconvex optimization problem that can be transformed into a nonlinear, convex 
optimization problem. 
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In our case, in the problem l |69l ), the constraints l !69b| > and l |69c| i contain functions that are non posynomial. Also, the 
variables in d69l > are not all positive, thus preventing a direct application of logarithmic transformation. In what follows, we 
first transform the problem l |69t into one equivalent in which the constraints involve functions that are all posynomial and the 
variables are all positive; and then we develop an algorithm for solving the equivalent problem. 

Let c = [ca, c;,, c, ] T s R 3 and 8 = [8 a , 8 b: 8 r ] T £ R 3 , such that ci > ^JPi/P and 8i = Pi + Ci for i = a,b, r. Note that the 
elements of 8 are all strictly positive. Also, for convenience, we define the following functions, for z = [z a , z b ] s Z 2 , 

Ki\' + \h b i\ 8 b c b ) + p(zl + zl)(\h al \ 2 {8 ) + |^fci| 2 (^6 + c fc)) 

( 2 2 \ / 2 2 \ 

\h a i\ z a S a c a + \h bi \ z b 8 b c b + h ai h M z a z b (S a c b + 8 b c a )J + N lz a + z b j 

V4(5, A 2 ,z) = A 2 (JV + P\h al \ 2 (8 2 a + cl) + P\h bt \ 2 {S 2 b + cl)) + 2P [z\ + z 2 b ) (\h al \ 2 8 a c a + \h bl \ 2 8 b c b ) 

+ P {\Ki\ 2 zi{8l + cl) + \h bt \ 2 z b (S b + cl) + 2h ai h bi z a z b (S a 8 b + c Q c b )) . (71) 

Let us now define the following functions, /i (6, A 2 ) = -01 (<*, A 2 , t), / 2 (6 , A 2 ) = (<*, A 2 , k), gi (<5, A 2 ) = V2 A 2 , t ) , 
g 2 (5,A 2 )=^(S,A 2 ,k), f 3 (S, A 2 ) = N + 2A 2 P\h rd \ 2 8 r c r , and g 3 (8, A 2 ) = A 2 (N + P\h rd \ 2 (8 2 r + c 2 )). 
It is now easy to see that the optimization problem $69\ can be stated in the following form. 



min A 2 (72a) 
<5, A 2 

91 {S,A 2 )-^ g 2 (8, A 2 ) g,{S,A 2 )-' ( ' Zb) 

~P~i [Pi 

— + Ci < Si < \ I — + Ci , i = a,b,r (72c) 

SeR 3 , eel 3 , A 2 el. (72d) 



The constraints d72b| l involve functions that consist of ratios of posynomials, i.e., are not posynomial — recall that a 
ratio of posynomials is in general non posynomial. Minimizing or upper bounding a ratio of posynomials belongs to a class 
of non-convex problems known as Complementary GP l29l . One can transform a Complementary GP problem into a GP 
problem using series of approximations. In order to get posynomials, we approximate the functions gi(8, A 2 ), g 2 (8, A 2 ) and 
33(5, A 2 ) with monomials, by using the following lemma ifTTl . 

Lemma 1: Let g{ S, A 2 ) = Y,j u j($i A 2 ) be a posynomial. Then 

^A 2 ),^,A 2 ) = n(^^) 7j . (73) 

Here, 7j = u 3 (8 (0) , A 2 0) )/g(8 (0) , A 2 0) ), Vj, for any fixed positive <5 (0) and A 2 0) then g(8 (0 \ A ( 2 0) ) = g(S (0) , A 2 0) ), and 
g(8^ , A 2 °^ ) is the best local monomial approximation to g(8^ , A 2 °^ ) near 8^ and A 2 °^ . 

Let gi(d,A 2 ), g 2 (8, A 2 ) and §3(8, A 2 ) be the monomial approximations of the functions gi(8, A 2 ), g 2 (8,A 2 ) and 
g 3 (8,A 2 ) obtained using LemmaQ] Using these monomial approximations, the ratios of posynomials involved in the con- 
straint fl72b| l can be upper bounded by posynomials. The optimal solution of the problem obtained using the convex approxi- 
mations is also optimal for the original problem d69l >, i.e., satisfies the Karush-Kuhn-Tucker (KKT) conditions of the original 
problem d69l >, if the applied approximations satisfy the following three properties 1301 , IfTTl : 



September 21, 2012 



DRAFT 



21 



1) gj{ S, A2) < §j(S, A2) for all S and A2 where g~j(8, A2) is the approximation of gj(8, A2). 

2) gj(S^°\ A^ ') = §j(S^ \ A^) where and Aj are the optimal solution of the approximated problem in the 
previous iteration. 

3) vgj(^°\ A^ ') = v5j(8^ \ A^ '), where vgj(-) stands for the gradient of function gj(-). 

Summarizing, applying the aforementioned transformations, we transformed the original optimization problem 469t first 
into a Complementary GP problem f72l and then into a GP problem by applying the convex approximations f73l . Finally, the 
obtained GP problem can be solved easily using, e.g., an interior point approach. More specifically, the problem of finding the 
appropriate preprocessing vector S for given integer vectors k and t can be solved using Algorithm A-2 hereinafter. 

Algorithm A-2 Power allocation policy for Rf°* as given by (|59]l 

1: Set 6 {0) to some initial value. Compute A 2 using S ( ' and set 1% = 1. 

2: Approximate g(8 (L2 \ A^' 2 ') with g(8 M , A<'' 2) ) around d (t2_1) and A^ 2_1) using {73j. 

3: Solve the resulting approximated GP problem using an interior point approach. Denote the found solutions as S^ 2 ^ and 
"2 ■ 

4: Increment the iteration index as 12 = 1^2 + 1 and go back to Step 2 using S and A2 of step 3. 
5: Terminate if ||<5 ( ' 2) - «5 (t ' 2 " 1) || < e or if A<'' 2) > A^ 2_1) . 



B. Compress-and-Forward at Relay and Compute at Destination 

The algorithms that we develop in this section to solve the optimization problem involved in the maximization of the sum 
rate given in Proposition[2]are essentially similar to those that we developed in the previous section. For brevity, we omit the 
details in this section. 

1) Problem Formulation: Recall the expression of as given by HT} in Proposition El The optimization 
problem can be stated as: 

(B) : max i min j log + I — — = ^ — 7 vf — I ' 

2 ^ \ snr||/3 s oH J a t -t|| 2 + (a t °a t ) F n d ) 

log + ( ... „ T ^r-)\, (74) 

\snr||/3 s oW a k - k\\ 2 + (a k o a k y n d J J 

where the distortion D is given by 

d _ iV 2 (l + snr||/3 s oh r || 2 ) iV 2 (snr(/3 s o h r ) T (/3 s o h d )f 
\h rd \ 2 P r \K d \ 2 P r (1 + snr||/3 s o h d \\ 2 ) 



and the maximization is over at, a k , (3 S such that < \/3 a \ < ^P a /P, < \/3 b \ < ^P b /P, and over the integer coefficients k 
and t such that |det(k, t) > 1. 

In order to compute R^m as given by l l74l . we develop the following iterative algorithm which optimizes the integer coef- 
ficients and the powers alternatively, and to which we refer to as "Algorithm B" in reference to the optimization problem 
(B). 
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Algorithm B Iterative algorithm for computing Rf°® as given by d74b 
1: Choose an initial feasible vector f3 s ^ and set t = 1. 

2: Solve J74b with /3 S = fi a lV> for the optimal k and t using Algorithm B-l and assign it to k^ 1 ' and t 
3: Solve J74b with k = k^ and t = t/'' for the optimal /3 a using Algorithm B-2 and assign it to f3 a . 
4: Increment the iteration index as 4 = i + 1 and go back to Step 2. 

5: Terminate if ||/3 s (l) - /V 1 ^ || < e, - - 1]| < e or if i&] < i& - !]• 



2) Integer Coefficients Optimization: Proceeding similarly as above, the problem of optimizing the integer 
vectors k and t for a fixed choice of the preprocessing vector (3 S can be written as 



mm 
k, t, e. 



s. t. 



Gi 



Gi > t T nt 

Gi > k T fik 



det(k,t) = |Mb - Ma | > 1 
k, t s Z 2 , 9i 6 R, 



(76a) 
(76b) 
(76c) 
(76d) 
(76e) 



where fl = (G T (GG T + N d ) _1 G - I 2 ) T (G T (GG T + N d ) _1 G - I 2 ) + ((GG T + N d ) _1 G) T N d ((GG T + N d ) _1 G). 
We reformulate the problem l l76l l into a MIQP problem with quadratic constraints [[2^ as before. We introduce the following 
quantities. Let a = [0,0,0,0, 1] T ; ai = a 2 = [0, 0, 0, 0, -1] T and a 3 = [0, 0, 0, 0, 0] T . Also, let b = [t a ,t b , k a , k b , Gi] T ; and 
the scalars c\ = c 2 = 0, and C3 = -1. We also introduce the following five-by-five matrices Fi, F 2 , and F3, where 





2n 
000 

The optimization problem d76b can now be reformulated equivalently as 



217 





2 



-2 




O000 



(77) 



mm 
b 



s. t. 



a b 

It t 

— b F^b + &i b < Cj i 

k s Z 2 , t £ Z 2 , Gi £ R 



1,2,3 



(78) 



It is easy to see that the matrices Fi and F 2 are positive semi-definite. However, the matrix F3 is indefinite, irrespective to 
the values of k and t. In order to transform the optimization problem ( 1781 into one that is MIQP-compatible (i.e., in which all 
the quadratic constraints are associated with semi-definite matrices), we replace the quadratic constraint l !76db with one that 
is linear, as performed in the previous section. 
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Finally, the MIQP optimization problem is given by, 



min 6i (79a) 

k,t,0! 

s. t. 6i > t T flt (79b) 

9i > k T flk (79c) 

-|«o^(l + k a +t b ) -K b T a (l + h + t a )\ < -1 (79d) 

k k 

— - 1 - ki < 0, — - + 1 + fcj < 0, i = a, 6 (79e) 

— - 1 - *i < 0, -— + l+ti<0, i = a,6 (79f) 
k, t £ Z 2 , k, t, k, t e R 2 ,9i £ R. (79g) 



The optimization problem l [79l can be solved iteratively using Algorithm B-l hereinafter. 



Algorithm B-l Integer coefficients selection for R^°° as given by (|74T > 
1: Initialization: set ij = 1, 

2: Use the branch-and-bound algorithm of 1271 . 1281 to solve for ©J , k' 11 ' and f/' 1 ^ with the constraint l !79dt substituted 

with -K a T b (l + k a + ti) + K b r a (l + k b + t a ) < -1. 
3: Update the values of k and t in a way to satisfy i65\ : and increment the iteration index as i\ = ti + 1. 
4: Terminate if exp(k (tl) ) « 1 + k (,l) and exp(t (tl) ) « 1 +t (tl) . Denote the found solution as 9" ,ln ' 2 . 
5: Redo steps 1 to 4 with in Step 2 the constraint d79d| > substituted with K a r b (l + k a + t b ) - n b r a (l + k b + t a ) < — 1. In Step 

4, denote the found solution as 9™ n ' 2 . 
6: Select the integer coefficients corresponding to the minimum among 9™"' 1 and 0™"' 2 . 



3 ) Power Allocation Policy: The problem of optimizing the power value (3 S for a fixed integer coefficients k, and 
t, can be written as, 



mm 

0,,oc t ,a k , © 2 



©2 



s. t. 9 2 > 



9 2 > 



snr||/3 s oH e*t-t||' + (at ° at) n d 



snr||/3 s o H T a k - k|| 2 + (a k o a fc ) T n rf 



D > 



N 2 (l + snr||/3 s o h r || 2 ) iV 2 (snr(/3 s o h r ) T (/3 s o h d )f 



|/i rd | 2 P r (l + snr||/3.oh d | 2 ) 



' P / P 

— < Pi < \J — , i = a,b 

f3 s eK 2 , 9 2 el. 



(80a) 

(80b) 
(80c) 

(80d) 

(80e) 
(80f) 
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Algorithm B-2 Power allocation policy for as given by d74t 

1: Set Sa^ ' ^ to some initial value and set L2 = 1 and 13 - 0. 
2: Compute G^ 2-1 '^, a t (t2_1 '' 3) and a ^ 2_1 ' t3) using Ss^ 2 ' 1 '^ . 

3: Approximate g(d (t2 '' 3) , 9^ 2 ' t3) ) with g(d s , (t2 ' t3) , 6^ 2 '' 3) ) around d^- 1 ' 13 ' and 9^ 2_1 '' 3) using {73. 

4: Solve the resulting approximated GP problem using an interior point approach. Denote the found solutions as Ss^ 2 ' 13 ^ and 
9(^3) 

5: Increment the iteration index as 12 = 12 + 1 and go back to Step 3 using 8 S and 02 of step 4. 
6: Terminate if ||<5 s (t2 ' t3) - <5 s (t2_1 ' t3) || < e or if 9 2 t2, ' 3) > Q[ L2 ~ 1>L3) and denote by S the final value. 
7: Increment the iteration index as L3 = 13 + 1, set 12 = 1, and <5 s ^ t2_1 '' 3 ^ = S and then go back to Step 2. 
8: Terminate if \R^[t 3 ] - -l]|Se or if iffli 3 ] < ifflis - 1]. 



As before, let c = [c a , q,] T s R 2 and S s = [<5 Q , <5fc] T e R 2 , such that > \fpjp and 5^ = ft + Ci for i = a,b. Note that the 
elements of <5 are all strictly positive. We can reformulate the optimization problem as, 



mm 
<5 S , OL t ,a k , e 2 



s. t. 



2 

fi(S,Q 2 ,a.t,a k ) 
gi(S,Q 2 ,at,a k ) 



< 1, 



g2(S,Q2,OL t ,a k ) 



< 1, 



/ 3 (a,e 2 ) 



< 1 



(81a) 
(81b) 



+ Ci < Si < 



S s s R 2 , C 6 R 2 , 9 2 6 



+ a , i = a,b 



(81c) 
(81d) 



The constraints <81bt correspond to the constraints d80bb . ( I80cb . and J80dt . These functions consist of ratios of posynomi- 
als, i.e., are not posynomial — recall that a ratio of posynomials is in general non posynomial. As before, we transform the 
Complementary GP problem into a GP problem using series of approximations. 

As there are several variables to optimize over simultaneously, the optimization is carried out in two steps. First, we 
optimize the power value S s using geometric programming with successive convex approximation in a way that is essentially 
similar to the above, for a fixed value of at and Next, we optimize the value of at and a k for a fixed value of <5 S . This 
process is repeated until convergence. More specifically, the problem of finding the appropriate preprocessing vector S s for a 
given integer vectors k and t can be solved using Algorithm B-2. 

V. Numerical Examples 

In this section, we provide some numerical examples. We measure the performance with the sum rate obtained under the 
constraint of symmetric rate, i.e., two times the maximum rate at which both users can communicate simultaneously with the 
destination, in a reliable fashion. We compare our coding strategies with those described in Section lTl-CI 

Throughout this section, we assume that the channel coefficients are modeled with independent and randomly generated 
variables, each generated according to a zero-mean Gaussian distribution whose variance is chosen according to the strength 
of the corresponding link. More specifically, the channel coefficient associated with the link from Transmitter A to the relay 
is modeled with a zero-mean Gaussian distribution with variance a 2 .; that from Transmitter B to the relay is modeled with 
a zero-mean Gaussian distribution with variance <r 2 r ; and that from the relay to the destination is modeled with a zero-mean 
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Gaussian distribution with variance a\ A . Similar assumptions and notations are used for the direct links from the transmitters 
to the destination. Furthermore, we assume that, at every time instant, all the nodes know, or can estimate with high accuracy, 
the values taken by the channel coefficients at that time, i.e., full channel state information (CSI). Also, we set P a = 20 dBW, 
P b = 20 dBW, P r = 20 dBW and P = 20 dBW. 




Fig. 2. Achievable sum rates under symmetric rate constraint. Numerical values are P = 20 dBW, ol = al = cr, 2 d = 20 dBW, 
and cr 2 d = cr b 2 d = dBW. 

Figure [2] depicts the evolution of the sum rate obtained using the so-called compute-and-forward at the relay approach, 
i.e., the sum rate R^m of proposition[T]as given by d22b : and the sum rate obtained using the so-called compress-and-forward 
at the relay and compute at the destination approach, i.e., the sum rate R^S of proposition[5]as given by fiTl i. as functions of 
the signal-to-noise ratio SNR = 101og(P/iV) (in decibels). Note that the curves correspond to numerical values of channel 
coefficients chosen such that the transmitters-to-relay links are of comparable quality to that of relay-destination link and 
are better than transmitters-to-destination links. For comparison reasons, the figure also shows the sum rates obtained using 
the classical strategies of Section Hl-CI i.e., the sum rate -R^i allowed by standard amplify-and-forward as given by dl4t . 
the sum rate RfJ m allowed by standard decode-and-forward as given by d!6l i. and the sum rate R^ m allowed by standard 
compress-and-forward as given by d21l i. All the sum rates are evaluated under the constraint of symmetric rate by the users. 

For the example shown in Figure[2] we observe that both strategies of proposition Q] and proposition [2] perform well. The 
strategy of proposition [2] achieves a sum rate that is larger than what is obtained using classical relaying techniques, for all 
ranges of SNR. Also, the strategy of proposition Q] outperforms standard AF and DF for all ranges of SNR. Figure [3] depicts 
the same curves for other combinations of channel coefficients, chosen such that the transmitters-to-relay links are of better 
quality than all other links, but the relay-to-destination link is not better than the transmitters-to-destination links. In this case, 



September 21, 2012 



DRAFT 



26 




Fig. 3. Achievable sum rates under symmetric rate constraint. Numerical values are P = 20 dBW, ai = al = 26 dBW, cr r 2 d = 14 
dBW, and <7 2 d = <r b 2 d = dBW. 



we observe that the curves show behaviors that are generally similar to those of Figure[5] Also, note that, the standard AF can 
provide sum rates that are better than those obtained using the strategy of proposition Q] at large SNR is precisely due to that 
the transmitters-to-destination links are of better quality than the relay-to-destination link for the example shown in Figure[3] 
(Recall that, although the strategy of proposition [T] does not require decoding the messages from the two users individually, 
its performance is limited by the bottleneck at the relay node nonetheless; and this explains why it can be outperformed by 
standard AF at large values of SNR as visible from Figure[3]i. Figure|4]shows similar curves for other combinations of channel 
coefficients, chosen such that the relay-to-destination link has better quality than all other links. 

Remark 2: Recall that, as we mentioned previously, the optimization Algorithm B associated with the strategy of propo- 
sition's non-convex. In the figures shown in this paper, the sum rate provided by this strategy, which is based on Wyner-Ziv 
compression at the relay and computing appropriate equations at the destination, are obtained by selecting only certain initial 
points for Algorithm B. For this reason, the sum rate offered by the coding strategy of proposition [2] i.e., R^°m, can possibly 
be larger than what is shown in the figures if one considers more initial points. Also, we note that even for initial points 
choices that yield a sum rate that is only as large as that obtained with standard compress-and-forward, there are advantages 
for using the coding scheme of proposition[2]mstead of regular CF. In particular, while regular CF utilizes the rather complex 
joint typicality decoding, or maximum-Likelihood decoding, the coding scheme of proposition[2]utilizes linear decoding, and 
thus is more easily feasible. 

Remark 3: The comparison of the coding strategies of proposition Q] and proposition |5] is insightful. Generally, none of 
the two coding schemes outperforms the other for all ranges of SNR, and which of the two coding schemes performs better 
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SNR 

Fig. 4. Achievable sum rates under symmetric rate constraint. Numerical values are P = 20 dBW, ai = al = 14 dBW, af d = 26 
dBW, and = <r b 2 d = dBW. 

depends on both the operating SNR and the relative strength of the links. For example, observe that while the strategy of 
proposition [5] outperforms that of proposition [Jin the examples shown in Figure S0 [3] and @] the situation is reversed for the 
example shown in Figure|5] Related to this aspect, recall the discussion in Remark[T] Let us denote by i?^m° D the achievable 
sum rate obtained by modifying the coding strategy of proposition [5] in a manner to get the first equation computed at the 
destination using only the signal received from the relay during the second transmission period, and the second equation 
computed at the destination using only the direct transmissions from the transmitters during the first transmission period. 
Figure[6] shows the improvement brought up by the aforementioned joint processing for the example depicted therein. 

We close this section with a brief discussion of the convergence speed of Algorithm A that we use to solve the optimization 
problem (A) given by i59i , as described in section lTV-AI Recall that the algorithm involves allocating the integer coefficients 
and the users' powers alternately, in an iterative manner. For a given set of powers, we find the best integer coefficients by 
solving a MIQP problem with quadratic constraints using the optimization software MOSEK. For a given set of integer- valued 
coefficients, we find the best powers at the sources and relay by solving a series of geometric programs by means of an interior 
point approach 1291 . 

In order to investigate the convergence speed of the proposed algorithm, we compare it with one in which the integer 
coefficients search is performed in an exhaustive manner and the power allocation is kept as in Section HV- A3 1 Note that, 
using this exhaustive-search algorithm, for the integer valued equations coefficients to be chosen optimally, the search can be 
restricted to the set of integer values that satisfy 1 1 k 1 1 2 < 1 + 1 1 h r 1 1 2 snr and 1 1 1 1 2 < 1 + 1 1 1 1 2 snr, since otherwise the allowed sum 
rate is zero 13 1 1 . Let Rfum denote the sum rate obtained by using the described exhaustive search-based algorithm. Figure [7] 
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Fig. 5. Achievable sum rates under symmetric rate constraint. Numerical values are P = 20 dBW, erj = of r = cr r j = 20 dBW, 
and al A = cr b 2 d = dBW. 



shows that the number of iterations required for Algorithm A to converge, i.e., yield the same sum-rate as the one obtained 
through exhaustive search, is no more than three. Also, we note that, in comparison, the exhaustive search-based algorithm is 
more largely time- and computationally resources consuming, especially at large values of SNR. Similar observations, that we 
omit here for brevity, also hold for Algorithm B. 

VI. Conclusion 

In this paper, we study a two-user half-duplex multiaccess relay channel. Based on Nazer-Gastpar compute-and-forward 
scheme, we develop and evaluate the performance of coding strategies that are of network coding spirit. In this framework, the 
destination does not decode the information messages directly from its output, but uses the latter to first recover two linearly 
independent integer-valued combinations that relate the transmitted symbols. We establish two coding schemes. In the first 
coding scheme, the two required linear combinations are computed in a distributive manner: one equation is computed at the 
relay and then forwarded to the destination, and the other is computed directly at the destination using the direct transmissions 
from the users. In the second coding scheme, the two required linear combinations are both computed locally at the destination, 
in a joint manner. In this coding scheme, accounting for the side information available at the destination through the direct 
links, the relay compresses what it gets from the users using Wyner-Ziv compression and conveys it to the destination. The 
destination then computes the desired two linear combinations, locally, using the recovered output at the relay, and what it 
gets from the direct transmission from the users. For both coding schemes, we discuss the design criteria and establish the 
associated computation rates and the allowed sum rate. Next, for each of the two coding schemes, we investigate the problem 
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SNR 

Fig. 6. Achievable sum rates under symmetric rate constraint. Numerical values are P = 20 dBW, ai = al = 14 dBW, af d = 26 
dBW, and = <r b 2 d = dBW. 

of allocating the powers and the integer-valued coefficients of the recovered equations in a way to maximize the offered sum 
rate. This problem is NP hard; and in this paper we propose an iterative solution to solve this problem, through a careful 
formulation and analysis. For a given set of powers, we transform the problem of finding the best integer coefficients into a 
mixed-integer quadratic programming problem with quadratic constraints. Also, for a given set of integer-valued coefficients, 
we transform the problem of finding the best powers at the sources and the relay into series of geometric programs. Comparing 
our coding schemes with classic relaying techniques, we show that the proposed strategies can offer substantial gains; and we 
illustrate this through some numerical examples. 
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